Low velocity quantum reflection of Bose-Einstein condensates 
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We studied quantum reflection of Bose-Einstein condensates at normal incidence on a square array 
of silicon pillars. For incident velocities of 2.5-26 mm/s observations agreed with theoretical predic- 
tions that the Casimir-Polder potential of a reduced density surface would reflect slow atoms with 
much higher probability. At low velocities (0.5-2.5 mm/s), we observed that the reflection probabil- 
ity saturated around 60% rather than increasing towards unity. We present a simple model which 
explains this reduced reflectivity as resulting from the combined effects of the Casimir-Polder plus 
mean field potential and predicts the observed saturation. Furthermore, at low incident velocities, 
the reflected condensates show collective excitations. 

PACS numbers: 34.50.Dy, 03.75.-b, 03.75.Kk 
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Quantum fluctuations of the electromagnetic field ex- 
ert forces on objects and are responsible for the attractive 
interactions between two neutral objects, e.g. an atom 
and a surface 0. Such interactions are typically weak 
and decay rapidly with increasing separation. Still, they 
are important for nanoscale devices 0, Q and when ul- 
tracold atoms are trapped close to a surface 0,01 El- O ne 
spectacular consequence of the Casimir-Polder potential 
is the prediction of total quantum reflection of very slow 
atoms from neutral surfaces: atoms incident on a sur- 
face at low velocity are accelerated toward the surface so 
abruptly that they reflect from the potential instead of 
being drawn into the surface UjBOlllJJ] ■ If high reflection 
probabilities could be realized, new atom-optical devices 
such as mirrors and cavities would be possible without 
the need of magnetic or optical fields. However, in a 
recent study of quantum reflection of Bose-Einstein con- 
densates, the reflection probability was limited to ~ 15% 
at low velocity 11] . A theoretical paper simulating quan- 
tum reflection of Bose-Einstein condensates could not ex- 
plain the low reflectivity 

In this work, we investigate the quantum reflection of 
Bose-Einstein condensates (BECs) from surfaces for ve- 
locities near and below 1 mm/s. To enhance quantum 
reflection, we use a pillared silicon surface, in the spirit of 
previous experiments with grazing-incidence neon atoms 
on ridged silicon 0,0,0. We observe quantum reflec- 
tion of Bose-Einstein condensates with probabilities of up 
to 67% for velocities of ~1 mm/s, corresponding to a col- 
lision energy of tjx 1.5 nK. We propose a simple model 
to explain how mean field interactions interfere with the 
reflection process and prevent the observation of higher 
reflection coefficients with BECs. Further, due to the 



greatly enhanced reflection coefficients, we observe col- 
lective excitations of the reflected condensate and inco- 
herent scattering between incident and reflected clouds. 

Bose-Einstein condensates of 23 Na atoms were pre- 
pared and transferred into a loosely confining gravito- 
magnetic trap, comprising a single coil and three exter- 
nal bias fields, as described in Ref. 0|. For typical load- 
ing parameters, condensates with N « 1 x 10 6 atoms 
were confined ~1 cm above the coil in a harmonic trap 
characterized by angular frequencies (w_i_, uj y , lo z ) = 2ir x 
(4.2, 5.0, 8.2) Hz, where directions (_L, y, z) are defined in 
Figure ^ At this point, lu± and u> y were adjuste d by 
changing the vertical bias field as described in Ref. jig . 
Typical densities in the trap were ~ 5 x 10 12 cm -3 and 
diameters were ~ 150 /«m. A silicon surface attached to 
a micrometric, motorized linear actuator was mounted 
~1 cm above the single coil. The position of the surface 
relative to the center of the coil was adjustable during 
the experiment as shown in Figure^. 

The surface used in this experiment, provided by the 
MIT Nanostructures laboratory, was a pillar structure 
etched into single-crystal silicon. The structure was cre- 
ated by Interference Lith ography (IL) and various sub- 
sequent etching steps [l7L fist llflj. Figured shows the 
final surface as an array of 1 /im tall, 50 nm diameter 
pillars spaced at 500 nm. Such a surface should provide 
a Casimir potential approximately 1% of the value for 
a solid Si surface. The quantum reflection efficiency de- 
pends inversely on the strength of the interaction, and a 
dilute surface is expected to exhibit enhanced reflection. 

Studying the reflection properties of the surface re- 
quires a controlled collision. After loading the condensate 
into the trap, the surface was moved to a desired distance 
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FIG. 1: Experimental schematic, (a) Atoms were confiried 
in a gravito-magnetic trap near a pillared Si surface. Atoms 
were accelerated towards the surface by displacing the trap- 
ping potential a distance d (greatly exaggerated) so that it 
was centered on the surface. The surface was mounted on 
a translation stage and could be removed at any point for 
imaging, (b) Scanning electron micrograph of the pillared Si 
surface used in this experiment. 



d from the trap center. By changing the bias field B± 
appropriately, a dipole oscillation centered on the sur- 
face was induced [llj. After waiting Tj_/4 = 27r/4u^ the 
atoms hit the surface with velocity v± = dto±. By vary- 
ing lu± between 27rx2 and 27rx4 Hz and d over 50 fim 
to 1 mm, velocities in the range of 0.5 to 26 mm/s could 
be studied. The reflection probability was calculated as 
the ratio of the average reflected atom number to the 
average incident atom number 26]. The reflection prob- 
ability, along with data previously collected for a solid 
silicon surface are shown in Figure El The pillared 
surface shows higher reflectivity over a wider range of in- 
cident velocity, as expected. The reflection maximum is 
67% for a velocity of 1.2 mm/s and reflection probabilities 
above 10% were measured at velocities up to 20 mm/s. 
Below ~3 mm/s, the reflection probability flattens near 
55%, qualitatively similar to the behavior of the solid 
surface where the reflectivity flattened near 12% in the 
same velocity range. 

We calculated a theoretical reflection probability of a 
single atom from the pillared surface using three numeri- 
cal simulations. The surface potentials of the Casimir- 
Polder form Ci/r A are obtained using Cf* = 6.2 x 
10~ 56 Jm 4 for bulk silicon [i(J and combining contribu- 
tions from both the pillar layer and the bulk substrate. 
Reflection probabilities were calculated by numeri cally 
solving the Schrodinger equation for a ID potential [2l|. 
We consider three averaging schemes for simulating the 
experiment: (1) we average the density of the material 
before calculating the potential, simulating the surface as 
a 1 /im thick overlayer of material with C4 = 0.01 xCf 1 
added to a semi-infinite slab of material with C4 =Cf*, 
(2) we calculate the 3D potential from the pillared struc- 
ture numerically by integrating over the regions of space 
containing material with C4 =Cf l [2^ and then average 
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FIG. 2: Reflection probability vs incident velocity. Data 
were collected in a magnetic trap with trap frequencies 
2tt x (2.0, 2.5, 8.2) Hz (squares) and 2tt x (4.2, 5.0, 8.2) Hz (tri- 
angles). For comparison, data from Ref. |lH| for reflection off 
a solid silicon surface are shown as circles. Incident and re- 
flected atom numbers were averaged over several shots. For 
clarity error bars for data below 5 mm/s are shown only on the 
inset plot, which has a different horizontal axis to emphasize 
the low velocity data. Systematic uncertainty in the veloc- 
ity due to residual motion is approximately 10%. Theoretical 
curves are described in the text. 



to obtain a ID potential, or (3) we assume that atoms 
follow a linear trajectory toward the surface (Eikonal 
approximation), and calculate the reflection probability 
from many points above the surface before averaging the 
reflection probability. The resulting reflection probabil- 
ity curves are shown in Figure |21 as numbered dashed 
gray lines. The predictions of model (1) and (2) are sim- 
ilar for the pillared surface. They show that the reflec- 
tion probability depends mainly on the the diluted pillar 
layer and only weakly on the bulk material underneath 
or the arrangement of the pillars. Model (3) should only 
be valid for high incident velocity, when the de Broglie 
wavelength X^b is significantly smaller than the surface 
structure. This is not the case in our experiment where 
XdB — l" m exceeds the spacing of the pillars. 

All calculations predict that the reflection probabil- 
ity approaches unity for low incident velocity. This is in 
contrast to our observation that the reflection probability 
saturates below 3 mm/s for both the pillared and solid 
surfaces. It was suggested that this saturation is due 
to low velocity excitations which smear out the conden- 
sate density. Although the reflectivity approaches unity, 
some reflected atoms would appear in a diffuse cloud 
which may fall below a detection threshold 121. How- 
ever, this could explain our previous results [ll| only 
when we assume a density threshold for detection of 
0.25 x no 10 12 cm~ 3 , where no is the central con- 
densate density, which is twenty times higher than the 
lowest densities we are able to detect [lij ]. 

There is a finite-size correction to the standard descrip- 
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FIG. 3: Mean field model for quantum reflection of conden- 
sates, (a) The trapped condensate provides a repulsive mean- 
field energy which is a constant away from the surface and, 
within the healing length f, drops to zero. The dashed curve 
shows this mean-field potential set to zero at infinity. This 
potential combined with the Casimir potential (dotted), cre- 
ates the composite potential (solid) which we use to model 
reflection in the presence of a condensate, (b) The reflection 
probabilities from the same the potentials for low velocities. 



tion of quantum reflection, but it is too small to account 
for our observations. For an incident atom cloud of size 
d, the smallest incident velocity is h/md, approximately 
0.2 mm/s for our parameters. We conclude that a single- 
particle description can not account for our low-velocity 
data and now discuss possible effects due to the conden- 
sate's mean field interaction. 

The mean field potential is taken to be that of a con- 
densate at rest with a fully reflecting wall as a boundary 
condition. The condensate's density decays towards zero 
at the wall over a characteristic length scale given by the 
healing length, £. The atoms at the edge of the conden- 
sate thereby acquire a velocity given by ~ h/m^ which 
is approximately equal to the speed of sound c. If the 
healing length is much larger than the relevant range of 
the Casimir-Polder potential, approximately 1 /xm as de- 
fined by the so called badlands region [lfj, one would 
assume that the mean field potential simply accelerates 
the atoms. Atoms leaving the condensate enter the re- 
gion of quantum reflection with an incident velocity ob- 
tained from mv /2 = U = mc 2 . This model would shift 
the single-atom quantum reflection curves by the velocity 
v = y/2c which is w 1.5 mm/s for our parameters. This 
shift is too small to explain the low reflectivity at our 
lowest velocities. Additionally, the assumption that the 
healing length is much larger than the distance at which 
quantum reflection happens is not valid for our data. 

In order to fully account for interaction effects, we cal- 
culate the quantum reflection probability using a com- 
posite potential which includes both the Casimir-Polder 



potential and the mean field potential (Fig^Ji). The re- 
sults now show a dramatic reduction of the reflectivity 
at low velocity, as shown in Figure^. At high velocities 
(> 3 mm/s), quantum reflection occurs close to the sur- 
face where the mean field potential plays no role. As a 
result, the predictions of the composite model are simi- 
lar to the single atom theory. As the velocity is reduced, 
the point of reflection moves outward, into the region 
where the mean-field potential "softens" the Casimir- 
Polder potential. At very low velocities (<0.1 mm/s), 
when the badlands region is far from the surface, the pre- 
dicted reflection resembles the reflection probability from 
the tail of the condensate rather than from the Casimir- 
Polder potential. We note that the model continues to 
predict perfect reflection at zero incident velocity, how- 
ever the approach is distinctly non-monotonic; the re- 
flection probability saturates and decreases precipitously 
near 1 mm/s before rising to unity. This model predicts 
well, without any free parameter, the velocities below 
which we have observed saturation of the reflectivity for 
both the solid and pillared surface as shown by the dot- 
dashed lines in Figure [21 Unfortunately, the data do not 
extend far enough into the low velocity regime to confirm 
the model's prediction of a sharp drop at low velocity or 
the ultimate asymptote to unity. 

The model does not include the effects of the mov- 
ing condensate, its observed collective excitations, or the 
distortion of the condensate wavefunction by surface at- 
traction or the loss of atoms to the surface. 

The calculated curves are not in quantitative agree- 
ment with the experimental data; the observed reflection 
probabilities are lower, even at high velocity. One pos- 
sible explanation is the further modification of the po- 
tential by stray electric fields, caused by sodium atoms 
deposited on the surface (adatoms). Recently, the par- 
tial ionization of rubidium adatoms by bulk silicon has 
been shown empirically to produce an electric field of 
several V/cm at 10 from the surface 0|. This elec- 
tric field, which falls off as 1/r 2 , will produce an ad- 
ditional potential, Vk(r) = — A/r 4 , which will reduce 
the reflection probability. To account for stray electric 
fields, we fit the high velocity data for the pillared (solid) 
surface using a potential Vt„t = —0.01 x Cf l /r 4 — A/r 4 
(Vtot = -Ct /r 4 -A/r A ). We find for the pillared (solid) 
surface a value of A of 0.02 x Cf l {Cf l ) corresponding 
to a stray field ~10 V/cm (~100 V/cm) at 1 fim for the 
pillared (solid) surface, smaller than the values measured 
in the rubidium experiment. If we combine the stronger 
surface potential with the mean-field potential we have a 
phenomenological model which is consistent with all our 
data, shown in Figure |3 as solid lines. It would be very 
interesting to test this model by varying the density over 
a large range and try to observe the predicted decrease of 
the saturation velocity for lower density. Unfortunately, 
we couldn't study quantum reflection at lower density 
due to rapid decrease of the signal-to-noise ratio. 
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FIG. 4: Menagerie of reflection effects. (a,b,c) As the incident 
velocity is reduced, the reflected condensate becomes increas- 
ingly excited, (d) By removing the surface at the moment 
of reflection, we can see both the incident (left) and reflected 
(right) condensates. The reflection probability is 30%. (e) 
The collision of the incident and reflected condensates pro- 
duces a strong s-wave scattering halo at low velocity, visible 
here Tx/4 after reflection. The surface is still present on the 
left in this image. Half of the halo is missing due to sur- 
face reflection or absorption, (f) With the surface removed, 
the scattered atoms remain in the trap after an additional 
half trap period, and appear reversed in position and veloc- 
ity. Field of view for images a, b, c is 540 /im and for d, e, f is 
800 (jm; the dashed line is the position of the surface (moved 
for imaging) at the moment of reflection. 



The higher reflection efficiency of the pillared surface in 
excess of 50% allowed us to study other aspects of quan- 
tum reflection of condensates at low velocity. Because 
of the finite size of the condensate, a standing wave is 
present in the condensate during the collision time which 
is inversely proportional to the incident velocity. As this 
time becomes comparable to the transverse and vertical 
trap periods, vortex rings, solitons, and other excitations 
may form, distorting the cloud |l2j . In our experiment, 
this velocity is approximately 2mm/s. At high velocities 
(> 4mm/s), we observe that the reflected condensate ap- 
pears, apart from diminished size and number, similar 
to the incident condensate. A condensed fraction and a 
thermal population, both present in the initial cloud are 
also present in the reflected cloud, shown in Figure^. As 
the incident velocity is reduced, as in Fig. 0Jd and c, the 
cloud develops a complex surface mode excitation |27|. 

Furthermore, we observe elastic s-wave scattering be- 
tween atoms in the incident and reflected condensates. 
S-wave scattering redistributes atoms in two colliding 
clouds with initial relative wavevector k evenly onto a 
sphere of radius \k\ [UEl; an i ma S e °f the atoms taken 
after a hold time T_l /4 will show the two clouds on oppo- 
site sides of the scattering sphere. In the present experi- 
ment, the reflected front part of the condensate collided 
with the still incident tail part (Fig.^Ji). The scattering 



halo was observed after sufficient hold time (Fig. _3e,f). 

We also performed the experiment using an aerogel 
surface. Aerogels are electrically insulating, randomly 
structured, silica foams with a density of ~2% of bulk 
silica |25) and should display reflection properties similar 
to the pillared surface. We were unable to observe quan- 
tum reflection above our detection threshold of ~2%, an 
effect we attribute to uncontrolled patch charges which 
strongly distort the Casimir-Polder potential and prevent 
efficient reflection. 

In light of the strong increase in reflection probability 
for the pillared structure, we want to discuss the ultimate 
limits of quantum reflection probability for condensates. 
We have found strong evidence that the presence of a 
condensate will distort the potential and prevent efficient 
reflection at low velocity. Our simple model predicts im- 
provements for longer healing lengths. Unfortunately, 
the corresponding reduction in condensate density would 
be a severe limit for atom optical devices based on quan- 
tum reflection. Another way to improve the reflectivity is 
by further reducing the density of the surface. Certainly 
it is possible to increase the pillar spacing. However, 
such widely spaced pillars will only dominate the poten- 
tial if their height is simultaneously increased, which is 
beyond the limit of current fabrication techniques. Sim- 
ilarly, narrower pillars may be possible, but not with the 
height required to dominate the potential. 
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